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A COMMUNICATION SYSTEM AND METHOD 
FOR PERFORMING FAST SYMBOL ESTIMATION 
FOR MULTIPLE ACCESS DISPERSE CHANNELS 



The present invention relates to wireless and wireline communication systems. More specif- 
ically, the present invention relates to an improved method for estimating symbol sequences 
of multiple sources sharing the same communication media. The invention is particularly 
useful for code-division-multiple-access (CDMA) receivers where the channel length is only 
limited to two neighboring symbols. 

In a scenario where m receivers are utilized, the received signal at each receiver is the sum 
of the output of channels driven by u different user signals plus the noise at the receiver. The 
channel between receiver i and user j can be denoted as hij. In general, ^ h qn if i ^ q or 
j 7^ n. Noise signals at different receivers are assumed to be independent of one another. If 
the user signals are actually CDMA signals and the largest duration of /i^-, i — 1, • • • , m and 
j = 1, • • , u is short enough compared to the symbol duration of those CDMA signals, the 
orthogonality of the spreading codes can guarantee the separation of CDMA signals at each 
receiver. Intersymbol interference (ISI) and multiple access interference (MAI) are negligible. 
When the largest duration of the signature waveforms is comparable with, or even larger 
than the symbol duration, ISI and MAI can become severe. The performance of detection 
made on symbols individually will not be good due to ISI and MAI. Such a situation calls 
for joint detection by which channel effects are modeled adequately and symbol sequences 
are estimated jointly. In minimum-mean-square-error (MMSE) based joint detection, which 
will be referred to as joint detection in the following, an inverse filtering process is used to 
remove intersymbol interference and multiple access interference. To do this, a coefficient 
matrix H formed from the channel responses is inverted. More precisely, the pseudo-inverse 
of this coefficient matrix: (H*H)~ l H* must be found. The main problem associated with 
joint detection is that the computation thus involved is intensive. The dimension of H is 
determined by the number of users, the length of spreading sequences, the length of symbol 
sequences, and the number of receivers involved (even with only one receiver, the effect of 
multiple receivers can be achieved by sampling a received continuous signal at a multiple 
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times of the symbol rate.). With a moderate number of users and a moderate number of 
symbols, the dimension of the coefficient matrix can be large. To speed up the inversion 
of the coefficient matrix, a Cholesky decomposition on H*H can be used. However, the 
computation amount involved in the Cholesky decomposition can still be substantial, as the 
requirement on matrix operations such as matrix inverses increases linearly with the number 
of symbols. It is therefore desirable to have a method to speed up the joint detection 
computation. 

Prior to this invention, a class of multiuser detectors have been developed. The most 
prominent ones among many others include R. Lupas and Verdu, "Linear multiuser detec- 
tors for synchronous CDMA channels", IEEE Trans, on Information Theory, 1(35):123-136, 
January 1989.; Z. Xie et al. "A family of sub-optimum detectors for coherent multiuser 
communications", IEEE Journal of Selected Areas in Communications, pages 683-690, May 
1990; and Z. Zvonar and D. Brady, "Suboptimum multiuser detector for syncrhonous CDMA 
frequency-selective Rayleigh fading channels", IEEE Trans, on Communications, Vol. 43 
No. 2/3/4, pp.154-157, February/March/ April, 1995, Angel M. Bravo, "Limited linear can- 
cellation of multiuser interference in DS/CDMA asynchronous systems", pp. 1435-1443, 
November, 1997. They either design various finite-impulse-response filters to process the re- 
ceived signals, which are in essence approximate solutions of the MMSE based joint detection 
solution, or use Cholesky decomposition to solve the equations involving the large coefficient 
matrix, e.g. Paul D. Alexander and Lars Rasmussen, "On the windowed Cholesky factoriza- 
tion of the time- varying asynchronous CDMA channel", IEEE Trans, on communications, 
vol. 46, no. 6, pp.735-737, June, 1998, which gives the exact MMSE joint detection result 
with high computational complexity, i.e. complexity of matrix operations such as matrix 
inversion increases linearly with the number of symbols. 

As in many communication systems, the length of ISI and MAI is often limited to sev- 
eral symbols, besides being Toeptliz or block Toeplitz matrix, the coefficient matrix is also 
banded. Fast algorithms developed for generic Toeplitz or block Toeplitz matrices, e.g. 
Georg Heinig and Karla Rost, "Algebraic Methods for Toeplitz-like Matrices and Opera- 
tors", Birkhauser, 1984, do not take advantage of the fact that the coefficient matrix is 
banded. The afore-mentioned Cholesky factorization based algorithm, while taking advan- 
tage of the banded nature of the coefficient matrix, does not take advantage of the fact 
that the coefficient matrix is also Toeplitz or block Toeplitz, which leads to considerable in- 
crease of computational demand as the number of matrix operations such as matrix inversion 
increases with the number of symbols to be estimated. Simpler versions of the method dis- 
closed in this invention, have been used to solve numerically Poisson equations, e.g. Roland 
A. Sweet, "A cyclic reduction algorithm for solving block tridiagonal systems of arbitrary 
dimension" , SIAM journal on Numerical Analysis, Vol. 14, pp. 706-720, September, 1977. 



SUMMARY OF THE INVENTION 

The present invention discloses a symbol detection system and method for multiple users 
in scenarios where the medium is shared by multiple users. The symbol detection system 
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and method involve formulating the received signals as the product of the inverse of a block 
triadiagonal matrix and a vector formed by the symbol sequences. The block tridiagonal 
matrix is derived from the channel responses. The present invention includes a fast algorithm 
to find the product. The computational complexity of the disclosed method is linear with 
the number of symbols and the requirement for matrix operations such as matrix inversion 
increases with the logarithm of the number of symbols. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a graphic illustration of a wireless communications setup for the disclosed technique. 

FIG. 2 is a graphic illustration of a wireline communications setup for the disclosed technique. 

FIG. 3 is an overview of the disclosed algorithm. 

FIG. 4 is a flow chart of the forward reduction procedure. 

FIG. 5 is a flow chart of the backward substitution procedure. 

While the invention is susceptible to various modifications and alternative forms, spe- 
cific embodiments thereof are shown by way of example in the drawings and will herein 
be described in detail. It should be understood, however, that the drawings and detailed 
description thereto are not intended to limit the invention to the particular form disclosed, 
but on the contrary, the intention is to cover all modifications, equivalents and alternatives 
falling within the spirit and scope of the present invention as defined by the appended claims. 

DETAILED DESCRIPTION OF THE DRAWINGS 
INTRODUCTION 

The basic problem under consideration is that of estimation of symbol sequences with inter- 
symbol interference and multiple access interference. To give a concrete setup of the joint 
detection problem under study, we consider a wireless CDMA system with multiple anten- 
nas (receivers) at the base station, and consider the joint detection problem of finding the 
symbol sequences from these users, i.e., the joint detection problem of the uplink signals 
received at the base station. As was made clear previously, by oversampling, the effect of 
multiple virtual receivers can be achieved with only one physical receiver, the applicability 
of the procedures which will be given is apparent for the joint detection problem of downlink 
signals also. In a code division multiple access system with J users, there are M (M > 1) 
antennas and hence M physical receivers at the base station, and the length of the spreading 
sequences is P. The received signal at each physical receiver is oversampled at M t times the 
chip rate thereby producing a plurality of received signals each sampled at the chip rate. 
In other words, for each physical receiver, the received signal is oversampled to produce a 
number of chip rate signals. Hence at the base station, the signal received at antenna i can 
be represented as follows: 
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n{n) = ]T hit{n) * c fc (n), i = 1, • • • M, 

where /^(n) is the signature waveform of the channel between the A-th user and the i- 
th antenna at the base station. The signature waveform can be identified either by blind 
methods or with training sequences, e.g. Tong et al, "Blind Identification and Equalization 
Based on Second-Order Statistics: A Time-Domain Approach" , IEEE Trans, on Information 
Theory, March, 1994. 

FIG. 1 shows one example of an embodiment where the invention may be used. A base 
station (102) receives signals from wireless terminals (104, 106, 108, 110). For some wire- 
less terminals(106, 108), there is more than one path between them and and the base station, 
as made clear by the reflective paths by reflectors (112, 114). Path length difference among 
multiple paths can destroy the orthogonality of signal waveforms of different transmitters 
which may have existed. Methods based matched filters and simple despreading are not able 
to offer satisfactory results under such circumstances. MMSE joint detection can offer much 
better results with heavier computation complexity. The disclosed method is able to reduce 
the computation amount requirement for implementing MMSE joint detection considerably. 

In the embodiment of FIG. 1, the symbol estimator of the invention may be comprised in 
the base station (102) for estimating symbols in uplink transmissions from wireless terminals 
(104,106,108,110). The symbol estimator of the invention may also be comprised in one 
or more of the wireless terminals (104,106,108.110) for estimating symbols in downlink 
transmission from base station (102). 

FIG. 2 shows another embodiment of the disclosed invention. A cable modem node 
(200) receives signals from multiple cable modems (206, 208, 210) which share the same 
cable medium (212) to communicate with the cable modem node (200). The channels from 
different cable modems to the cable modem node (200) will be different. In the embodiment 
of FIG. 2, the symbol estimator of the invention may be comprised in the cable modem node 
(200). Similar to the situation shown in FIG. 1, the disclosed technique can also be applied 
to speed up the joint detection method that combats the dispersive channel effects. It is not 
difficult to see that the invented technique can be applied in a reversed way, i.e., the cable 
modems (206,208, 210) may include the symbol estimator of the invention and may receive 
the signal from the cable modem node (200). 

FIG. 3 shows one embodiment of the method used by this invention to estimate symbol 
sequences of multiple users. The base-band signals are collected and sampled from multiple 
receivers (302). In (304), At each receiver, the channel response between a transmitter 
and that receiver is used to filter the sampled received signal at that receiver. As there are 
J channel responses between all the transmitters and that receiver, there are a total of J 
filtered outputs produced from the sampled received signal at receiver. All the outputs are 
multiplexed to form a single data vector, according the transmitter indices, in an ascending 
order. Next all the multiplexed data vectors from all the receivers are added up to form one 
single vector Y\ (318). 

Channel responses h m j(n), m = 1, M, j — 1, J (312) between transmitters 
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and receivers are identified either by training sequences or other methods which are not 
specified in this invention. Space-time correlation matrices T;, i = d, d - 1, * • • , — (d — 
1), -rf, are constructed from those channel responses (312). Channel description matrices 
Ai, Bi, and Dj are block Toeplitz matrices, the blocks of which are defined by Tj, i = 
^ d — l 7 . . . 7 —(d - l), — d 7 which are constructed to form the description matrices (316). 
Forward reduction steps (306) are initialized with the second set of matrices and vectors, i.e. 
Ax, B x ,Di and Y x . In the course of forward reduction, using the bisection approach, some 
intermediate matrices which will be used in the backward substitution steps in (308) are 
stored. At the beginning of backward substitution steps (308), the intermediate solution X 3 
is found. Solved unknowns and intermediate matrices which have been stored in the forward 
reduction steps are used to find other unknowns. The output of the backward substitution 
steps (308) is demultiplexed (310) to produce the estimate of symbol sequences sent by 
transmitters. 

FIG. 4 is a flow chart showing the forward reduction procedure. 
FIG.5 is a flow chart showing the backward substitution procedure. 



DATA FORMATION 

In this disclosure, cases where signature waveforms of arbitrary duration are considered. In 
the following, we will write the received signals as the convolution of signature waveforms 
and oversampled symbol sequences in a structural way, which facilitates the derivation of 
the fast algorithm. If the signature waveform duration is limited to d + 1 symbol duration, 
then we can pad zeros to h iik (n) so that the length of /^(n) is precisely (d + 1)P. Assume 
that hmj(n)j 1 < m < M, 1 < j < J are column vectors, which are identified with either 
training sequences or other means, which is not specified in this invention (312). We first 
define some notation. We divide /w,-(n) into (d + 1) parts: 



where /^ mj '\ j = 1, d + 1, are P x 1 vectors. We put together the signature waveform 
vectors: 



h m i 



h 



771 J 



4' 



(ml) 



(ml) 



(mJ) 



(mJ) 



■(ml) 
l d+l 



(mJ) 
*d+l 



H (») 1 



H 



(m) 
d+1 
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where 



■(ml) 



H£ m) is a P x J matrix, k = 1, • • • , d + 1. 

Put the received signal at receiver m into a vector fashion: 



L + d< 



f r«(l) 
r-(2) 

[r m ((I + d)P) J 



(*») 



H 



(m) 



Hi' 



Tj(™) 



X 2 (l) 

xdL) 
x 2 (L) 

. x i( L ) . 



X 



By stacking the received signals: Rm, m = 1, • • • , M, one over another, we have 



r * i 




"H x * 




' Wi ' 






H 2 
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. Rm . 
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The objective of joint detection is to find X. If 

{L + l)PM>LJ, 

then normally H is a full rank matrix. 
The zero forcing solution is given by 



X = (H*H) _1 H*i2 = (£ H^H™)" 1 £ H m i^. 

m=l m=l 

The MMSE (minimum mean square error) solution is given by 

M M 

X = (H*H + a 2 I)- l H*R = ( £ H^H m + ^I)" 1 £ H^, 



m=l 



m=:l 



where a 2 is the variance of the noise. We define a matrix: 



M 



T := H*H + v!=Y, H mH m + «I, 



m=l 



and a vector 



M 



m=l 



where v can be either 0 or a 2 depending on whether the zero-forcing solution or the MMSE 
solution is in need. H^i^, m = 1, • • • , M are formed in (304). Y is formed in (318). In 
both zero-forcing and MMSE solutions, the computation of X can be separated in two steps: 

1. Find Y. 

2. Find T~ l Y. 

In either case, T is a block banded matrix, defined in the following way, 

T 0 T_ d 



T_ fl 



(5) 



0 iLJxLJ 



where 



T.= EE( H SUrH£U +i + ^K T-^TJ, » = 0, -1 -d. 

m=l Jb=0 

And T»,i = d, d - 1,« • -(d - 1), -d are called space-time correlation matrices (314), 
delta(-) is the Kronecker function. 

We can represent T as a block tridiagonal matrix, if d divides L: 



1. the diagonal block A is given by 

A = 

T*-i 



(d-i) 



2. the supra-diagonal block B is given by 



B = 



T_ 



d JdJxdJ 



We can write 



as 



TX = Y 



L < 



A B 
B* A B 



B* A B 
B* A 



X = Y. 



(6) 



(7) 



FAST JOINT DETECTOR 



DERIVATION 

The system of equations as given in (7) can be treated as the special case of a more general 
formulation defined by a triplet [A, B, D]: 



M Bi 
B*i A< Bi 



B\ Di 



Xi — Yi, 



(8) 



where A,, Bj, and D f are J x J matrices, which are called channel description matrices 
(316), 



Xi 



(9) 
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(10) 



with Xij and Y^ being J x 1 vectors. 

We linearly combine equations for the following cases: 



1. for 2k = 2 (i.e. = 1), 
2* - 1 : -B^A" 1 



2A;: 
2A; + 1 



I 

-B.A- 1 



x 
x 

X 



A; Bi 

B*i Ai Bi 

B*; ^ Bi 



^i,2fc-l 
-X"t,2fc+2 J 



*i,2A-l 
^i,2A: 



■A 



(11) 



and we have 

[0 -B*i A t -1 Bi + Ai - BiA _1 B*i 0 -B^B, ] 
= -B*iAf 1 y r i, 2 fc_i + Yi» - BiAT^+i, 



or 



[ -B*iA- x Bi + Ai - BiA~ x B*i -Bi A^B* ] 
= -B*iA~ 1 li, 2 jk_i + Y iak - BiA- 1 ^^!. 



-Xi,2Jfe+2j 



-Xi,2fc+2 



(12) 



(13) 



2. For every even index 2&(4 <2k<L { - 3), we can linearly combine the 2k - 1 st , 2k th , 
and 2k + 1 st equations in system of equations (8): 



2k - 1 : -B'iA* 1 
2k: I 



2fc + l 



— BiA 



-l 



x 
x 

X 



B*i Ai Bi 

B*i Ai Bi 

B*i Ai Bi 



^i,2*-2 
^i,2Jb+2 



\ 



Yi,2k-l 
Yi t 2k 
Yi,2k+1 . 



J 

(14) 



and we have 



[ -B* i A t " 1 B*i 0 -B^A-'Bi + Ai-BiA-^i 0 -BiA^B* ] 
= — B*iA" 1 3 / i j 2jt_i + Vi )2 jfc - &iAj' l Yi£k+u 



^z,2Jb+l 
-Xi,2*+2 



(15) 
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or 



[ -B* i A l " 1 B% -B*iAr l B< + Ai - B<A- l B*i -B^A^B* ] 



= -B*jA i 1 li,2fc-l + *i,2A — BjAi x Yi y 



^i,2Jk-2 
-X"i,2A 
. -Xi,2fc+2 



,2Jb+l- 



(16) 



3. Depending on whether Li is odd or even, we need to study two different cases. 



(a) If Li is odd, the last even index 2k between 1 and Li is Li — 1. Let Li = 2j + 1, 
so fc = j . We combine the last three equations 



2j-l 
2j: 



-B*iAr x x 
I x 



2j + l: -BiD: 



we have 



-l 



B\ Ai Bi 

B*i Ai Bi 
B\ Di 



Yi j2 j-l 
*i,2j+l 



\ 



(17) 



[ -B* i Af 1 B* i 0 -B*iA~ l Bi + Ai — BjD~ 1 B* i o] 



-X"i,2j-1 
-Xi,2j 
-X»,2j+1 J 



(18) 



or 



[ -B^A-^^ -B^A^B; + Ai - BiD- l B*i ] 
— B* i A~ l Yi i 2j-i + Yi^j — B*iDj 1 li j 2j+i* 



X h2j 



(19) 



(b) If is even, the last even index 2k between 1 and Li is Li, Let Li = 2j. We 
combine the last two equations: 



2 j - 1 : -B'iAT 1 
2j: I 

we have 



x 
x 



B% Ai Bi 
B*i Di 



Xi,2j-2 
-Xi,2j-1 



Y i,2j 



, (20) 



[-B%AT l B% 0 -B^A-'Bi + Di ] 



Xi$j-2 
Xi,2j-l 



= — B*jA i 1 i^,2j-i + ^i,2i? 

(21) 
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or 



* -B* i AT 1 B* < — B*iA^ 1 B» + Di ] 



^-B\A-% 2j ^ + Y iy2j . (22) 



To facilitate a structural description of the newly-generated equations, we can define the 
following: 



If Li is odd, we define 
and 

If Li is even, we define 
and 

We also define 



A i+1 — B^A^Bj + Ai - B^Af 1 B* £ , 
Bi + i := — B^A"^. 

D i+1 := A, - BiD^B'i - B^A^B*, 
Li+i :— -{Li — 1). 
D i+1 := D.-BV^B,, 



Li+i :- -Li. 



(23) 

(24) 
(25) 
(26) 
(27) 
(28) 



Xi+i,k : — Xi^ki 1 < 2k < Li. 

Aj + i> Bj + i, and D i+i are the new channel description matrices. 

Substitute the above newly-defined variables into equations (16) and (13) for 1 < 2k < 
Li — 2, and equation (19) if Li is odd, or equation (22) if Li is even. We obtain the following 
system of equations: 



B*j+i Ai + i B» + i 



B*i+i Aj + i Bj + i 
B*j+i Di + i 



(29) 



We call the system of equations given in (29) as the i + 1 st system of equations. Hence a 
system of equations of reduced dimension involving only Xij of even indices as unknowns is 
obtained in the cancellation process, and Yi+\ is obtained. It can be seen that the coefficient 
matrix defined by the triplet [A i+1 , B^+i, D i+ i] is still a block-tridiagonal matrix, even 
though of smaller dimension than the coefficient matrix defined by by the triplet [A*, Bj, D,]. 
Hence the same procedure can be again applied to this new system of equations. With the 
execution of the procedure each time, the size of the system of equations is reduced by at 
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least a half. If we start from the 1 st system of equations as given in (10), by executing 
this procedure not more than log2{L\) times, we can generate the 2 nd , 3 rd , • • • systems of 
equations. Eventually, at some s, the 5 - 1 th system of equations would be given either as 



B%_i A 5 _i B s _i 
B* s -i D s _x 



Xs-l — Ys 



or 



A 5 _i B s _i 
B%_i D 5 »i 



(30) 



(31) 



(30) and (31) can be reduced to the s** 1 system of equations involving one uxl subvector: 
X $i i as unknown, respectively through (17) and (20): 



we have 



(32) 



(33) 



The just- found X Si x(ov X 5 _i,2 as In the s — 1 st system of equations) can be substituted into 
the s — l 5 * system of equations. In general, for the i th system of equations, from solution 
of the i + 1 st system of equations, X iy j of even indices are already known, and Xij of odd 
indices are still unknown. We have three cases to consider. 



1. Consider the first equation in the I th system of equations: 



Xi v 2 



X 



(34) 



where X^X^i) has been found from the solution of the % + 1 st system of equations. 
We have 

X^^A-Hni-B^u). (35) 



2. Consider the 2k + 1 st equation in the i th system of equations, = 1, 2, etc.: 



2* + l: [B% A* B^ 



Xi 



(36) 



where the values of X^^i+i,*) and Xi )2 jt+2(^i+i,&+i) have been found from the solu- 
tion of the i + 1 st system of equations. We have 



or 



Xi^k+i = A~ l (Yi£k+i — &*iXi+i,k — BjXi+i^+i)- 



(37) 
(38) 
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3. If Li is odd(Ii = 2k + 1), then X i>Ls is found from the last block equation: 



[B% D £ ] 



X i,Li-l 
X i,li 



(39) 



with 



Xi^t = D 4 (Yi, L . - B*iX i)L .-i). 



(40) 



Consequently for the i th system of equations, we just solve the block equations corre- 
sponding to odd indices. Hence all the unknowns in the i th system of equations are found. 
The above procedure can be recursively applied to the i - I s *, i - 2 nd , • • • systems of equa- 
tions until the unknowns in the system of equations as given in (7) are all found. 

Apparently, the main computation is on the right-side of the equations. The update of 
the triplets involves light computation. 

In summary, the solution of the block tridiagonal system of equations as given in (10) is 
executed with two procedures. In the first procedure, the size of the system of equations is 
reduced, which is called forward reduction (306), a flow chart of which is shown in FIG. 4. 
In the second procedure, solved unknowns are substituted into larger systems of equations 
and more unknowns are solved, which is called backward substitution (308), a flow chart of 
which is shown in FIG. 5. 



COMPUTATIONAL PROCEDURE 

The whole procedure is naturally separated into two parts: the forward reduction procedure 
and the backward substitution procedure. 
The forward reduction procedure: 

1. Ai := A, B 2 := B, D l ~ A, Y x := Y, L x := L, i := 1. 

2. Find the inverse of A*: A" 1 . 

3. Qi := BiA~ l and E; := B*^ 1 . 

4. A i+l := A* - EfBi - (E^Bi)*- 

5. B i+1 := — QjBj. 

6. • If Li is odd, 

(a) P^BA'S 

(b) D m :=A i -P i B%-E i B i , 

(c) L M := ^ 

(d) Yi+ lJt := -EiYut-x + Y ij2k - Q-Yw+u l<k< 

(e) Y mM+l : = -E,y iiLi _ 2 + Y iM . x - P^. 
• If Li is even, 
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(a) D i+1 ^Di-EiBi, 

(b) X i+1 := f , 



(c) li+i,* := -EiYi, 2k -i +Y it2k - Qi^jt+i, l<k< L i+l , 

(d) *i+i,L<+i := -EiY i>Li -i + Y itLi . 

7. Store A^, D^\, Qj, E t , Pj (if it is calculated), and Fj. 

8. i := i + 1. 

9. Repeat 2 to 8 until Lj = 1. 
10. $:=i. 

A flow chart of the above-stated procedure is given in FIG. 4. 
The backward substitution procedure: 

1. Solve D S X S = Y 3 . 

2. i := s - 1. 

3. -Xi,2fc := -Xi+i^) 1 < & < Li + i. 



6- 7 := j + 2, if j < I/,- go to 5. 

7. i := i - 1, if * > 0 go to 3. 

8. Xi is then the estimate of X. 

A flow chart of the above-stated procedure is given in FIG. 5. 

Although the system and method of the present invention has been described in connec- 
tion with the preferred embodiment, it is not intended to be limited to the specific form set 
forth herein, but on the contrary, it is intended to cover such alternatives, modifications, 
and equivalents, as can be reasonably included within the spirit and scope of the invention 
as defined by the appended claims. 




if j = 1. 

if 1< j < Li. 

if j = 
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CLAIMS 



What is claimed is: 

1. A method for estimating at least a subset of symbol sequences from received signals 
at multiple receivers in a multiple receiver communication system wherein each of said 
multiple receivers receives a sum of signals from a plurality of transmitters sharing 
the same communication medium, wherein said signals are formed by a sequence of 
contiguous symbols, each of which occupies a fixed time duration, and samples said 
received signals with a time interval that is the inverse of an integer multiple of said 
fixed time duration, wherein said integer multiple is no smaller than the number of said 
signals divided by the number of said receivers, wherein the channels of said signals 
have finite impulse response and are different from one another, said symbol estimation 
method comprises 

(a) sampling said received signals at multiple receivers with a time interval equal to 
a j~ fraction of said fixed time duration, wherein F, an integer is no smaller than 

wherein J is the number of said transmitters, M is the number of said 
receivers, M t is a positive integer and generating Q sampled received signals from 
each of said received signals wherein Q is an integer no larger than M t \ 

(b) constructing data vectors from sampled received signals by filtering them with 
channel responses between said transmitters and said receivers and multiplexing 
the outputs; 

(c) obtaining the maximum length of the channel responses of users at multiple re- 
ceivers, which is denoted as d+ 1, and constructing a channel response matrix for 
each receiver, 

(d) constructing space-time correlation matrices with said channel response matrices; 

(e) constructing a second set of matrices and vectors computed from said space-time 
correlation matrices and said data vector; 

(f ) applying a recursive bisection approach on said second set of matrices and vec- 
tors to generate an output sequence, wherein said recursive bisection approach 
comprises 

forward division steps which recursively compute a first series of matrices and 
vectors initially based on said second set of matrices and vectors and yield vectors 
in said first series of vectors and matrices with dJ[Li/2\ rows where dJLi is the 
number of rows of the last element in said first series of vectors and matrices, 
until a pre-defined stopping criterion is satisfied, where [xj means the smallest 
integer which is not smaller than x, 

backward substitution steps which first construct an intermediate solution with 
the last of the first series of matrices and vectors and recursively compute, in a 
reverse order of said forward division steps, a second series of vectors with said 
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first set of matrices and vectors until exhausting all of said first set of matrices 
and vectors; 

(g) demultiplexing the vector from the last in said second series of vectors; whereby 
said symbols of said signals are found from said output sequence. 

2. The symbol estimation method in claim 1 wherein said channels are convolutions of 
physical channels and their corresponding pulse shaping filters of said channels. 

3. The symbol estimation method in claim 1 wherein said channels are convolutions of 
physical channels and their corresponding spreading codes of said signals in a code- 
division-multiple-access communication system. 

4. The symbol estimation method in claim 1 wherein said channels are convolutions of 
physical channels, their corresponding spreading codes of said signals, and their cor- 
responding pulse shaping filters in a code-division-multiple-access communication sys- 
tem. 

5. The symbol estimation method in claim 1 wherein said channel matrix comprises chan- 
nel filter responses. 

6. The symbol estimation method in claim 1 wherein the vectors in said second set of 
matrices and vectors are constructed by adding up said data vectors. 

7. The symbol estimation method in claim 1 wherein each of said channel response ma- 
trices is constructed for each said receiver by putting the channel responses between 
said transmitters and said receiver as columns in a matrix according to the user index, 
zeros are padded to said channel response matrices if necessary so that the number 
of rows in all said channel response matrices is the smallest integer multiple of of the 
maximum channel response length. 

8. The symbol estimation method in claim 1 wherein said space-time correlation matrices 
are constructed by adding up the correlation results of channel response matrices. 

9. The symbol estimation method in claim 1 wherein the matrices in said second set of 
matrices and vectors comprises three matrices A x , Bi, and D\ which are constructed 
by first finding the maximum length of channel filter responses of all users, padding 
zeros to channel filter responses whereby all of them are of d + 1 symbol durations 
wherein d is the least integer possible, constructing d time correlation matrices for 
each receiver with d matrices obtained from dividing the channel responses at that 
receiver into d + 1 parts, constructing d spatial-time correlation matrices for multiple 
receivers by adding up the time correlation matrices for all receivers with the same 
time shift, constructing A L according to a block-Toeplitz structure with spatial-time 
correlation matrices 0, • • • , d- 1, constructing Bi according to another block-Toeplitz 
structure with spatial-time correlation matrices -1, • • • , -<2, and finally assigning Ai 
to Di. 
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10. The symbol estimation method in claim 1 wherein said first series of matrices and 
vectors comprises a series of three dJ xdJ channel description matrices which are A* 
the ordinary diagonal block matrix, the superdiagonal block matrix, D* the special 
diagonal block matrix, and a right-hand vector which is divided into d J x 1 blocks 
which are indexed from one; 

11. The symbol estimation method in claim 1 wherein said forward reduction step further 
comprises putting said second set of matrices and vectors as the first of said first series 
of matrices vectors and computing recursively new channel description matrices from 
the latest channel description matrices and a new right-hand vector with blocks from 
linear combinations of blocks of even index and their neighboring blocks in the last 
right-hand vector. 

12. The symbol estimation method in claim 1 wherein said pre-defined criterion is that the 
length of the latest generated right-hand vector is dJ or that the superdiagonal block 
matrix in the latest of said first series of matrices and vectors becomes a zero matrix. 

13. The symbol estimation method in claim 1 wherein said backward substitution step 
further comprises solving the equations as defined by the latest of said first series 
of matrices and vectors, putting latest channel description matrices and the solution 
as the first of said second series of matrices and vectors, assigning the blocks of the 
vector of the latest in said second series of matrices and vectors to the even blocks of 
a new vector and computing the odd blocks of the new vector from these blocks and 
corresponding right-hand vectors. 

14. The symbol estimation method in claim 1 wherein said d spatial-time correlation ma- 
trices Ti, i = 0, * • • , —d, are constructed by 

M d+i 



m=l n=l 



where the channel response matrix Hi at receiver m, m — 1, • • • , M, consisting of 
(d + I) J rows, where d+ 1 is said maximum channel length, and J columns, which 
are the channel responses from transmitters to the receiver m, is related to H^, n = 
1, •••,<*+ 1 by 

where <$(•) is the Kronecker function, v is the variance of noise for minimum-mean- 
square-error joint detection and 0 for zero-forcing joint detection. 

15. The symbol estimation method in claim 1 wherein said second set of matrices and vec- 
tors consists of channel description matrices A 1? Bi, &i and Y x where A x is composed 
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of d x d blocks, and the (i, j) block is T_j +i if j > i, and Tiy +i if j < i; B x is composed 
of d x d blocks, and the (i, j) block is T, d+J _j if i > j, and Q JXJ if i < j; D x is A£; 

and Y\is % Ri where J2j is said signal obtained from receiver i. 

16. The symbol estimation method in claim 1 wherein said second set of matrices and 
vectors are B^A' 1 , B*Af \ B<Df 1 if is odd, X i+X , and Yj, for backward substitution 
step z, z = 5, 5-1, L 

17. The symbol estimation method in claim 1 wherein said forward division step comprises 
for forward division step i, i = 1, • • • , s, 



= -B^A^B* + Ai - B i A~ 1 B* i 

= -B^A-^i 

= -B^A^B** 

= ^A-^.x + ^-B.A" 1 ^^, = L^J 



(41) 



D i+ i 
Li+i 



if ^ is odd; 



if Li is even. 



Li+i 



= Ai - BiDf x B*i - B* i A l " 1 B i 
= 1(^-1) 

:= Dj - B*j Af 1 B i 
— It- 



(42) 



(43) 



18. The symbol estimation method in claim 1 wherein said pre-defined criterion is Li+i = 1 
or B i+ i — 0, 5 is set to i + 1 and the forward reduction steps stop. 

19. The symbol estimation method in claim 1 wherein intermediate solution step comprises 
solving B $ X S = Y $ . 

20. The symbol estimation method in claim 1 wherein said backward substitution step i, 
i = s, • • • , 1, comprises 



- - A i ^i-Yj+u, 

*i f 2*+i - ^yiakn-K^iX^-^BiX^x^ l<2k + KL h (44) 
X iA = D~% Li - D^B^i+x^ if L t is odd. 
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21 . A method for estimating at least a subset of symbol sequences from a 
received signal at a receiver in a communication system, the method comprising: 

(a) oversampling said received signal, thereby producing a plurality of 
received signal sample sets; 

(b) constructing data vectors from the received signal sample sets 

(c) obtaining the maximum length of the channel responses of users at 
multiple virtual receivers, which is denoted as d+1, and constructing a channel response 
matrix for each receiver, 

(d) constructing space-time correlation matrices with said channel response 
matrices; 

(e) constructing a second set of matrices and vectors computed from said 
space-time correlation matrices and said data vector; 

(f) generating an output sequence by applying a recursive bisection approach 
on said second set of matrices and vectors, wherein said recursive bisection approach 
comprises: 

forward division steps which recursively compute a first series of matrices and 
vectors initially based on said second set of matrices and vectors and yield vectors in said 
first series of vectors and matrices with dJ[Li/2] rows where dJLi is the number of rows 
of the last element in said first series of vectors and matrices, until a pre-defined stopping 
criterion is satisfied, where [x] means the smallest integer which is not smaller than x, 

backward substitution steps which first construct an intermediate solution with the 
last of the first series of matrices and vectors and recursively compute, in a reverse order 
of said forward division steps, a second series of vectors with said first set of matrices 
and vectors until exhausting all of said first set of matrices and vectors; 

(g) detecting symbols in the received signal by demultiplexing the vector from 
the last in said second series of vectors; whereby said symbols in the received signal are 
found from said output sequence. 



22. The method of claim 21, wherein each of said multiple receivers receives 
a sum of signals from a plurality of transmitters sharing the same communication 
medium, wherein said signals are formed by a sequence of contiguous symbols, each of 
which occupies a fixed time duration, and samples said received signals with a time 
interval that is the inverse of an integer multiple of said signals divided by the number of 
said receivers, wherein the channels of said signals have finite impulse response and are 
different from one another. 
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23 . The method of claim 2 1 , wherein said o versampling said received signal is 
performed at multiple receivers with a time interval equal to a 1/P fraction of said fixed 
time duration, wherein P is an integer no smaller than J/MM t , wherein J is the number of 
said transmitters, M is the number of said receivers, M t is a positive integer and 
generating Q sampled received signals from each of said received signals wherein Q is an 
integer no larger than M t . 



24. The method of claim 21 , wherein said constructing includes constructing 
the data vectors from the sampled received signals by filtering them with channel 
responses between said transmitters and said receivers and multiplexing the outputs. 
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Abstract 

A method for estimating symbol sequences from multiple transmitters which communi- 
cate with common receivers sharing the same medium is disclosed. This invention can be 
used in wireless communications and wireline communications such as in cable modem sys- 
tems and digital subscriber loops to speed up joint detection of symbol sequences of multiple 
transmitters to achieve low error probability. For convenience, only applications of this in- 
vention to wireless communications employing multiple antennas are described herein. The 
requirement for matrix operations like matrix inverses increase with the logarithm of the 
number of symbols instead of the number of symbols as with techniques based on Cholesky 
decomposition. 
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FORWARD REDUCTION PROCEDURE 
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